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ABSTRACT 

We present a study of the star-forming properties of a stellar mass-selected sample 
of galaxies in the GOODS NICMOS Survey (GNS), based on deep Hubble Space 
Telescope imaging of the GOODS North and South fields. Using a stellar mass selected 
sample, combined with HST/ACS and Spitzer data to measure both UV and infrared 
derived star formation rates, we investigate the star forming properties of a complete 
sample of ^1300 galaxies down to log = 9.5 at redshifts 1.5 < z < 3. Eight percent 
of the sample is made up of massive galaxies with M* ^ IO^Mq. We derive optical 
colours, dust extinctions, and ultraviolet and infrared star formation rates (SFR) to 
determine how the star formation rate changes as a function of both stellar mass and 
time. Our results show that SFR increases at higher stellar mass such that massive 
galaxies nearly double their stellar mass from star formation alone over the redshift 
range studied, but the average value of SFR for a given stellar mass remains constant 
over this ~ 2 Gyr period. Furthermore, we find no strong evolution in the SFR for our 
sample as a function of mass over our redshift range of interest, in particular we do not 
find a decline in the SFR among massive galaxies, as is seen at z < 1. The most massive 
galaxies in our sample (log M* ^ 11) have high average star formation rates with 
values, SFRuv,corr = 103 ± 75 M Q yr -1 , yet exhibit red rest- frame (U — B) colours at 
all redshifts. Wc conclude that the majority of these red high-redshift massive galaxies 
are red due to dust extinction. We find that A2800 increases with stellar mass, and 
show that between 45% and 85% of massive galaxies harbour dusty star formation. 
These results show that even just a few Gyr after the first galaxies appear, there are 
strong relations between the global physical properties of galaxies, driven by stellar 
mass or another underlying feature of galaxies strongly related to the stellar mass. 
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1 INTRODUCTION 

One of the least understood aspects of galaxy formation and 
evolution is the triggering and regulation of star formation in 
galaxies, and how this varies across galaxy mass and time. 
Galaxies are formed in star formation events, and under- 
standing how and when star formation occurs is critical for 
understanding the evolution of galaxies. One intriguing as- 
pect of the star formation in galaxies is that the star for- 
mation history appears to vary strongly as a function of 
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galaxy mass, at least at z < 1.5. Even within the local uni- 
verse we observe that the most massive galaxies are qui- 
escent ellipticals, while lower mass galaxies are undergoing 
significa nt star formation in either irregular or spiral galax- 
ies (e.g. IConselicd 120061. One key way to address galaxy 
formation and evolution directly is to understand how the 
nearby galaxy population was put into place and evolved 
from higher redshift galaxies which we can now observe in 
nearly complete mass-selected samples up to 2 = 3 (e.g. 
iDaddi et al1l2007l : IConselice et al.ll201l| y 

Observational results at redshifts z < 1.4 reveal that 
the most massive galaxies (M» ^ 10 n M o ) appear to 
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be dominated by old stellar populations with negligible 
amounts of star formation, and are largely as s embled by z ~ 
1 llMcCarthv et al] 12004: iDaddi et all [20051: ISaracco etaD 
l2005l ; lLabbe et alll2006l ; IConselice et alj|2007h . These galax- 
ies therefore must have formed very quickly inside of matter 
over-densities, within a few billion years after the Big Bang, 
then stopped forming stars, presumably as cold gas sup- 
plies diminished. This is in broad agreement with analyses 
of stellar popula tions seen in nearby massive ellipticals (e.g. 
IZhu et all l2010h , a nd is consistent with the broad tre nds 
seen in colours (e.g. lFaber et al.ll2007l ; iKriek et alj|2008h . 

Following up on this, several studies have shown 
that star formation rates (SFRs) correlate with stellar 
mass l|Daddi et al.ll2007l ; iNoeske et a"fll2007t iPannella et all 
2009) up to z — 2 and that the most massive galax- 
ies appear to shut off their star formation at ear- 



lier times than low-mass galaxies (e.g 


. Cowie et all 1996; 


Brinchmann & Ellis 2000; Bauer et all 12005; Feulner et al. 


2005; Perez-Gonzalez et all 120051; 1 Juneau et all 


2005; 


Bundv et all 20061; Paoovich et al. |200( 


; Zheng et al. 


2007; 


Schiminovich et all 120071: iDamen et al 


120091; iBauer et al. 


2010f). This effect of star formation ending in higher mass 



galaxies before lower mass systems is known as "downsiz- 
ing." The origin of downsizing is still not understood and 
must be the result of the removal of cold gas from a galaxy 
or the heating of this gas by some feedback process. Possi- 
bilities for the origin of this feedback include: active galactic 
nuclei (AGN), supernova, energy from a merger, or a com- 
bination of these. Standard galaxy formation models have 
difficulty reproducing downsizing as well as the properties 
of red der galaxies observed at high redshifts (e.g. I Cole et all 
2000), and properties of modern massive red galaxies that 
harbour old stellar populations. 

Furthermore, we know very little about how star for- 
mation and galaxy mass assembly rates vary as a function 
of galaxy mass and time at high redshift. While there is evi- 
dence that galaxy mergers vary with stellar mass at high red- 
shifts, such that the most mass ive galaxies are undergoing 
the most merging at z > 1 fe.g. lConselice et al.ll20ol . l2008l ; 
iBluck et all 12003 ). it is largely unknown how other modes 
of formation, such as cold gas accretion or in-situ star for- 
mation, are driving the formation of galaxies at early times. 
A deep, and at least moderately large area, survey complete 
in stellar mass is needed to address this point. 

There exists some evidence for how massive galaxies 
at high redshifts form stars, and the nature of their stel- 
lar population content, based on their optical and near in- 
frared colours and limited spectroscopy at 2 ~ 2. Mas- 
sive galaxies in the distant universe have been observed to 
be often very red, creating what is possibly the origin of 
the "red sequence" f o und in nearby ma ssive galaxies (e.g. 
IConselice et all 120071 ; IKriek et all 120081 ). Usually it is as- 
sumed that these red galaxies seen in the distant universe 
consist of old stellar populations. However, it is possible that 
these massive galaxies are red due to dust, or a combination 
of older stellar populations and dusty star formation. It is 
therefore possible that the dust content of massive galaxies 
in the universe has been underestimated, and a full consen- 
sus of the dust properties of high redshift massive galaxies 
is not yet known. 

There is evidence that dusty star formation is com- 
mon in producing red colours in lower redshift galaxies at 



z < 1. For example, by u sing a Jf-selected sample of galax- 
ies, ICimatti et all |2002l ) found that among K s < 19.2 Ex- 
tremely Red Objects (EROs) at z ~ 1, 50% are old pas- 
sively evolving galaxies while the other half are dusty star- 
forming gal a xies. In a more recent study of EROs at z ~ 1, 
iKong et all l|2009h confirm these results, showing a passive 
fraction of ~48% at K < 19.2, in an otherwise active star- 
forming population. 

There is also evidence that more distant galaxies are in 
fact redder due to dusty star format ion, rather than from 
old stellar populations. For example, IPannella et al.l (|2009l ) 
examine a sample of star-forming BzK galaxies from the 
COSMOS survey at z > 1.5 and derive dust attenuation 
factors from both UV and radio wavelengths. They find a 
strong agreement between the two methods of measuring 
dust content, and also show that for their sample of galaxies 
with median redshift z ~ 1.7, dust extinction increases with 
stellar mass. However, what is missing from these studies 
is a stellar-mass-selected sample of galaxies to systemati- 
cally investigate these trends, and to determine how the star 
formation rate is changing with stellar mass and redshift. 
Furthermore, we are also interested in further investigating 
the role of dusty star formation in the observed trend that 
massive galaxies are red at high redshift, while lo wer mass 
systems tend to be bluer (|Griitzbauch et al]|2010l) . 

One way to address these issues is to directly probe 
the star formation rate evolution of high redshift galaxies 
as a function of stellar mass and time. This allows us to 
determine directly the assembly rate for these systems at 
early times, and how this mode of assembly varies as a 
function of the already-established relationshi p with stel- 
lar mass. This has been done to some extent l|Daddi et all 
120071 ; INoeske et al1l2007l ). but has been largely limited to 
examining the most massive galaxies, and furthermore only 
a subset of these which are star forming. 

In this paper, we study the properties of a mass- 
selected, complete sample of galaxies from th e GOODS 
NICMOS Survey (GNS) i|Conselice et all 1201 ll ). focusing 
the analysis on star-forming properties of galaxies with 
log M* ^ 9.5 in the redshift range of 1.5 < 2 < 3.0. In 
§ [5] we describe the observations and sample selection. We 
describe our methods for determining physical properties of 
the galaxies in § [21 including a discussion of the reliability of 
SFR indicators in this redshift range. Section [3] discusses 
the trends of star-forming properties with colour, stellar 
mass, and redshift. We summarise our conclusions in § [SJ 
Throughout this work we adopt an Qm = 0.3, Qa = 0.7, 
Ho — 70 km s _1 Mpc _1 cosmology, use a Salpeter initial 
mass function, and provide all magnitudes in the AB sys- 
tem unless otherwise noted. 



2 DATA AND SAMPLE 

In this section we introduce the data we use from the 
GOODS-NICMOS Survey (GNS) and describe the deriva- 
tion of photometric redshifts and stellar masses. We then 
discuss the completeness of the survey at the redshifts of 
the study and finally introduce the supplementary infrared 
data. 
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2.1 The GOODS-NICMOS Survey 

Our galax y sample comes from the GOODS-NICMOS Sur- 
vey (GNS: IConselice et alJliuTll ). The GNS is a large Hub- 
ble Space Telescope program conducted with the NICMOS 
NIC3 instrument, to obtain deep //-band and J-band imag- 
ing within the GOODS North and South fields (GOODS: 
Grea t Observatories Origins Deep Survey; iDickinson et al.l 
l2003h . Using 180 HST orbits, a total of 60 pointings were 
observed, centered on the most massive galaxies (M* ^ 
10 n M Q ) at 1.7 < 2 < 2.9. The massive galaxies on which 
the pointings are centred were selected based on previous 
colour-selected samples within the GOODS fields that were 
designed to find distant galaxies, including BzA'-selected 
galaxies, Extremely Red Objects (EROs), and Distant Red 
Galaxies (DRGs). While the primary field selection was done 
in terms of the massive galaxy selection through these three 
primary colour criteria, to optimise the field placement, cat- 
alogues of Lyman-break BM/BX objects as well as high red- 
shift drop-outs and sub-mm galaxies, were also used. For a 
full de scription of the sample selection see IConselice et al.l 
<|201ll ). 

The main observations for the GNS survey were car- 
ried out using the NICMOS-3 camera on HST to obtain 
3 orbits per pointing for each image. With this exposure 
time, we achieved a F160W (//-band) depth of 26.8 mag 
(AB) at 5a using a 0.7 arcsec-diameter aperture. The field 
of view of each image is 51.2 arcsec on a side with a pixel 
scale of 0.1 arcsec pixel -1 . The FWHM for our images is 
rough ly 0.3 arcs e c. The data reduction is described in detail 
by iMagee et~afl (l2007f) who detect sources in the //-band 
using SExtractor (Bcrt in fe Arnouts! ll996). The full GNS 
catalog is derived to include every non-stellar object identi- 
fied in the //-band imaging over the total area of the survey: 
45 arcmin 2 . Overall, we find a total of 8298 objects in the 
GNS area. 

We supplement the near-infrared observations by 
matching sourc es to the GOODS/ ACS optical catalog 
(|Giava lisco 2004j). Photometry in B, V, i and z is available for 
sources down to a limiting magnitude of B ~ 28.2. After ac- 
counting for the known 0.3 arcsec offset in declination in the 
ACS v2 catalog, we perform source matching between the 
NICMOS///-band and ACS/z 8 50-band catalogs. The mean 
separation between optical and //-band coordinates their 
ACS counterpart is ~ 0.07 ± 0.15" for both the North and 
South fields. For the original matching between samples, a 
radius of up to 2" was allowed, but in the resulting sample, 
only 18 of the matches have a separation greater than 0.5". 
We do not find a strong correlation between catastrophic 
failures in photometric redshift and having a matching ra- 
dius separation greater than 1 arcsec, or 1.5 arcsec separa- 
tion so we have left these objects in the sample. 



2.2 Redshifts 

Where possible, we use spectroscopic redshifts published 
in the literature for GNS galaxies, otherwise we use our 
own measured photometric redshifts. Spectroscopic red- 
shifts of sources in the GOODS-N field were compiled by 
iBarger et ail <j2008T >. whereas the GOODS-S field spectro- 
scopic redshifts are taken from the FIREWORKS compi- 
lation l|Wuvts et al.l [20081 ). We matched these catalogs to 



the full GNS photometric catalog with a matching radius 
of 2", finding very few multiple matches, and a mean sep- 
aration between photometric and spectroscopic sources of 
0.11 ± 0.06" in the GOODS-N field and 0.13 ± 0.05" in 
the GOODS-S field. In the full GNS sample, there are 537 
spectroscopic redshifts for sources in GOODS-N and 369 in 
GOODS-S. 

For the remaining GNS galaxies, photometric red- 
shifts were obtained by fitting template spectra to the 
BVizH photometric data points, following the standard y 2 
minim ization procedure usi ng HYPERZ ([B olzonclla et al.l 
2000) as described fu lly in iGriitzbauch et al.l (|2010h and 



Conselice et al.l (|201 lh . The synt hetic spectra used by HY - 
PERZ are constructed with the iBruzual fc Chariot! (|2003t ) 
evolutionary code representing roughly the different mor- 
phological types of galaxies found in the local universe. We 
use five template spectra corresponding to the spectral types 
of E, Sa, Sc and Im as w ell as a single burst s cenario and a 
reddening law taken from lCalzetti et all (|200QT ) . We compute 
the most likely redshift solution and corresponding probabil- 
ity in the parameter space of age, metallicity and reddening. 

With the filters we use to determine photometric red- 
shifts the reddest rest-frame wavelength we probe is /3-band 
at z ~ 3. We cannot use longer wavelength observations to 
determine photometric redshifts because there are no near- 
IR observations redder than //-band that have the same fi- 
delity and depth as the BVizH bands we use in this study. 
The ground based /('-band observations available for this 
field are not nearly as deep as our NICMOS //-band obser- 
vations. While we have IRAC data for our sources, we do not 
use these data to determine photometric redshifts or stellar 
masses due to issues with the PSF and contamination from 
neighbouring galaxies. Furthermore, the rest-frame B-band 
gives us a good anchor f or m easuri ng stellar masses, as is 
shown bvlBell fc de Jonl (e.g. l200ll ) at lower redshifts, and 
iBundv et alj (120061 ) for higher redshifts. 

As with any wide-field study of galaxies at z > 1, the 
accuracy of photometric redshifts presents a challenge when 
interpreting observational results and cannot be neglected. 
We define the reliability of photometric redshift measures 
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z P hot)/{l + z apec ). We compare 



the average error ((Az/(1 + z))) and rms scatter (a) as well 
as the fraction of catastrophic outliers, i.e., galaxies with 
|Az/(l + 2)| > 0.3. Using HYPERZ, we find that sources 
with a high probability (P > 95%) have (Az/(1 + z)) = 
0.03, with a scatter of a = 0.045 (356 out of 906 galaxies 
with P > 95%). Considering all sources regardless of their 
probability we find similar values: (Az/(1 + 2)) = 0.01 and 
a = 0.06. We therefore consider all photometric redshifts 
regardless of their probability, which significantly increases 
our sample size, but only slightly increases the photometric 
redshift uncertainty. The fraction of catastrophic outliers for 
the full sample is ~ 6%. 

If we consider only galaxies in the redshift range of 
1.5 < 2 < 3 (as we do in this study; see Section 12. 4p . 
we find that 133 of the 1282 members of the final cata- 
log have spectroscopic redshifts. We find an average offset 
(Az/(l + z)) = 0.06 and arms of uaz = 0.10, with a fraction 
of catastrophic outliers of 20%. 

The most reliable photometric redshifts in the 2 > 1 
regime have so far come from studie s of galaxie s that 
have extensive multi-band coverage (jllbert et al.l 2009; 



4 A. E. Bauer et al. 



ICardamone et all l2010h . For example, lllbert et al.l (|2009l ) 
derive photometric redshifts for the COSMOS sample of 
galaxies out to z ~ 1.25 using 30 unique filters and account- 
ing for the emission line contribution to broad-band magni- 
tudes. They find an accuracy of (Az / (1 + z)} = 0.007. Using 
t he same technique f or the zCOSMOS galaxies at 1.5 < z < 
3. lllbert et ail (|2009h find an accuracy of (Az/(l + z)) = 0.06 
and a 20% fraction of catastrophic failures. This is similar 
to the range found in this study, but note that they define 
catastrophic outliers as objects with acc uracies worse than 
15% ( as opposed to the 30% used here). ICardamone et all 
(|2010i ) use 18 medium-band filters for a sample of galax- 
ies between 1.2 < z < 3.7. They determine an accuracy of 
photometric redshifts of (Az/(1 + z)) — 0.02 and a rms of 
o-az = 0.03. 

We caution in the paper that the reliability of our pho- 
tometric errors is not as high as with studies that are able 
to use 18 to > 30 filters to determine photometric redshifts. 
However we note that when comparing our photometric red- 
shift accuracies with a study that use a similar number of 
filters to derive photometric redshifts, we find that the reli- 
ability is in good agreement. Using the CDFS. lDamen et al.l 
(2009) find a median (Az/{l+z)) = 0.03 for their full sample 
and (Az/(1 + z)} = 0.079, for galaxies with z > 1. Again we 
stress that the unique aspect of this study is ability to mea- 
sure stellar masses at a deeper level than previously given 
our deep H-band HST data. Redshifts are a critical aspect 
of this study however, and we revisit this issue in detailed in 
Section 14.21 where we discuss the possible effects of the pho- 
tometric redshift errors on all our results through several 
simulations. 

2.3 Stellar Masses and Restframe Colour 

Stellar masses and restframe colours are determined from 
multicolour stellar population fitting techniques using the 
same catalog of five broad band data points used to deter- 
mine photometric redshifts for all GNS galaxies. A detailed 
description of how stellar masses a nd rest-frame (U — B ) 
colo urs are derived can be f ound in IConselice et al.l (|201lT ) 
and lGriitzbauch et al.l (|201ot ). and is summarised in the fol- 
lowing. 

To calculate these masses and colours we construct a 
grid of model spe c tral e nergy distributions (SEDs) from 
iBruzual &: Charlotl (|200St ) stellar population synthesis mod- 
els, assuming a Salpeter initial mass function and a varying 
star formation history, age, metallicity, and dust extinction. 
The star formation history is characterized by an exponen- 
tially declining model of the form tf>(t) = ip a exp(— t/r) with 
the star formation timescale ranging from r = 0.01 to 10 
Gyr and the age of the onset of star formation ranging from 
to 10 Gyr. The dust content is parametrized by the V- 
band optical depth with values of ry — 0.0, 0.5, 1, 2, and the 
meta llicity ranges from 0.0001 to 0.05 l|Bruzual fc Charlotl 
l2003h . 

The magnitudes resulting from the model SEDs are 
then fit to the observed photometric data of each galaxy 
using a Bayesian approach. For each galaxy, a likelihood 
distribution is computed for stellar mass, age, and abso- 
lute magnitude for each possible star formation history. We 
chose to compute rest-frame (U — B) colours, since the wave- 
length range of the U- and B-bands are covered best by the 
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Figure 1. [z — H) vs. (B — z) diagram for GNS galaxies with 
log M* > 9.5 and 1.5 < z < 3.0. Large circles show galaxies with 
log Af* ^ 11 and crosses identify galaxies also detected at 24 
jim. The solid diagonal line and dashed horizon tal line represent 
the constraints of criteria of lDaddi et al. for star-forming 

galaxies. The other lines represent the evolution of single stellar 
populations of given ages as observed at 1.5 < z < 3.0. Galaxies 
with ACS/-B-band magnitudes fainter than 28.2 (AB) are given 
this value as a limit and shown in grey with arrows. 

observed optical and //-bands. The peak of the likelihood 
distribution is adopted as the galaxy's stellar mass and Mu 
and Mb absolute magnitudes, while the uncertainty of these 
values is determined by the width of the distribution. While 
parameters such as age, e-folding time, dust extinction, and 
metallicity are not accurately fit due to various degeneracies, 
the stellar masses and colours are robust. From the width 
of the probability distributions we determine typical errors 
for our stellar masses of 0.2 dex. There are additional uncer- 
tainties from the choice of the IMF and due to photometric 
errors, resulting in a total random error of our stellar masses 
of ~0.3 dex, roughly a factor of two. The rest-frame (U — B) 
colours have errors of ~ 0.2. 



2.4 Sample Selection and Completeness 

For this work, we use a complete sample of galaxies selected 
by stellar mass. In this way we can perform an analysis of 
the assembly of all galaxies down to log M* = 9.5 at high 
redshift without obvious biases towards star-forming galax- 
ies. A detailed analysi s of stellar mass functi on for the GNS 
Survey is presented in lMortlock et all l|201ll ). Here we sum- 
marise the method of selection by stellar mass. 

To determine the limiting stellar mass for the GNS, 
we compute the expected completeness limits from the 5 a 
magnitude limit of the survey and the mass-to-light ratios 
(M/L) of simple stellar populations (SSP). The M/L ratio 
at the rest-frame wavelength corresponding to the observed 
NICMOS///-band at z ~ 3 is determined for a maximally 
old galaxy using the models of IWorthevI l| 19941 ). The M/L 
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ratio is then multiplied with the luminosity corresponding 
to a z = 3 galaxy at the detection limit of Hab = 26.8 to 
obtain the limiting stellar mass. In this way, we calculate a 
stellar mass complete ness limit of log M* = 9.2 at z = 2.75 
iMortlock etail (120111 ). 

The deep limiting //-band magnitude easily exceeds 
that of the expected upper bound for dusty bright sub- 
millim eter galaxies of 23.3 magnitudes found bv lFraver et al.l 
so we are confident that we are not missing the dusti- 
est galaxies due to attenuation. But even so, to test if we 
are biased against low mass dusty star-forming objects we 
examine with stellar population models if we would find, for 
example, a very low mass dusty galaxy at z = 2.5. If this 
object is composed of 50% old stars and 50% 100 Myr stellar 
population, this object would be 2.82 magnitudes brighter in 
M_g than a galaxy composed of a purely old stellar popula- 
tion formed at z — 10. Thus, for this object to be undetected 
at the very faint end of our completeness would require a 
dust extinction of 3.16 magnitudes, which is an amount of 
dust extinction exceeding what is commonly found in the 
lowest mass galaxies in our sample. That is, we would be 
detecting these extremely dusty low mass galaxies at our 
depth, but are not likely because they do not exist in any 
great numbers. 

From the original GNS catalog of 8298 galaxies, we se- 
lect a sample of galaxies complete down to log M* = 9.5 
and 1.5 < z < 3.0, which corresponds to a sample size of 
1282 galaxies within a survey volume of 2.3 x 10 5 Mpc 3 . We 
find that 8% of the galaxies are massive (log M* ^ 11) and 
10% have spectroscopic redshifts. Our final sample is the 
largest HST-based dataset with rest-frame optical imaging 
of mass-selected galaxies complete down to log M* = 9.5 
over 1.5 < z < 3. 

To gain a preliminary understanding of the properties of 
this mass-selected sample, we look at t he optical an d near - 
infrared two-colour plot developed by iDaddi et al.l (120041 ) 
as a technique to determine whether galaxies at z ~ 2 
show properties of star-forming or passive systems (the BzK 
method). With the ass umption that (H — K) = 0.25 on av- 
erage at our redshifts l|Brammer fc van Dokkum|[2007h . we 
show the corresponding BzH colour-colour space in Fig- 
ure [T] for our complete sample. Black points show galaxies 
with 5 a detections at all three magnitudes. The full mass- 
complete sample was selected from H-band observations to 
H ^ 26.8. The zsso-band depth is z = 27.5 which results 
in a 98.5% (1263/1282) detection rate among the full mass- 
selected sample. The 5 a detection limit in the /?435-band is 
B = 28.2 which detects 87.5% (1120/1282) of the galaxies 
from the mass-selected sample. 

Non-detections are assigned the limiting value (B = 
28.2) and are shown as grey arrows in Figure U Small points 
show all GNS galaxies with log M* ^ 9.5 while large cir- 
cles represent massive galaxies, with M„ ^ 1O 11 M0. Crosses 
represent galaxies also detected at 24 |im with /24 > 30/iJy. 
The solid diagonal black line and dashed horizontal black 
line show the BzH regions defined by (z — H) — (B — z) < 
-0.2 - (H - K) AB and (z - H) > 2.5 - (H - K) AB . The 
variou s tracks show passively evolving iBruzual fc Chariot! 
(2003) stellar population models as they would be observed 
between z = 1.5 — 3 from a single burst occurring at the 
ages indicated at the upper left. 

Most of the GNS galaxies in Figure [Tj sit to the upper- 



left of the solid diagonal line, indicating th ey are star- 
formin g galaxies according to the criteria of IDaddi et al.l 
(2004). More massive galaxies tend to have redder (z — H) 
and (B — z) colours, with almost all the massive galaxies 
showing (z — H) > 2. Galaxies detected at 24 fim are spread 
throughout the colour-colour space. In contrast, there are 
only two galaxies (detected in all three bands) within the 
triangular region to the upper right, indicative of passively 
evolving stellar populations. It is notable that neither of 
these galaxies are detected at 24 fim. Among the galaxies 
with limiting B-band magnitudes (grey) in the passive re- 
gion, only one is detected at 24 /j,m. The dot and dash-dot 
evolutionary tracks in the passive region of Figure [T] show 
galaxies of age 900 Myr and 1.4 Gyr, respectively, as they 
would appear between z = 3 and z — 1.5 when the universe 
was between 2.2 and 4.4 Gyr old. This shows that using 
colour-colour diagrams does not work well for selecting red 
galaxies since you need extremely deep B-band observations. 

2.5 Supplementary Infrared Data 

The Spitzer GOODS Legacy program l|Dickinson et al.l 
2003) provides deep infrared observations covering the 
GOODS fields, which we use to calculate infrared proper- 
ties for GNS galaxies. Since the resolution of the Spitzer 
Multiband Imaging Photometer (MIPS) 24 fim images is 
6" FWHM, we use the Spitzer Infrared Array Camera 
(IRAC) data to perform source matching within a radius 
of 1". IRAC data were reduced with the pipeline version 
S18.7.0 and then mosaicked with M0PEX, a package devel- 
oped at the Spitzer Science C enter for image proces si ng. A s 
described in Appendix A in IPerez- Gonzalez et"afl ( 2008), 
IRAC sources detected with SExtractor (down to a 5-er 
level of /(3.6/im) ~ 24.0 mag) were cross-correlated with 
UV-to-NIR images and photometry was measured in consis- 
tent apertures to build complete SEDs. A complete list of 
the extended data set s used , and their properties, is given in 
IPerez-Gonzalez et al.l (|2008l ). For the work presented in this 
paper, the redshif ts were fixed to the values used throughout 
this paper and in lGriitzbauch et al.l (|2010l ). 

The photometric method also includes a deblending al- 
gorithm for the IRAC and MIPS images which takes the 
known positions of the blended sources obtained from opti- 
cal ACS images and the point spread functions (PSFs) for 
the different images, an d obtains separated fluxes for the 
blended sources (see also iGrazian et all j|2006D ). 

This paper makes use of the MIPS observations of GNS 
galaxies to calculate infrared luminosities and IR star for- 
mation rates. The MIPS 24 /im data are those obtained as 
part of the GOODS Guaranteed Time observations and the 
Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), 
of both North and South fields. The final mosaic has a 
5-a detection level of 30 /xJy within the GOODS area. 
C atalogs were made wit h a PSF-fitti ng method described 
in IPerez-Gonzalez et al.l (|2005l . 120081 ). The MIPS catalog 
has an 80% completeness limiting flux of 30 /iJy for the 
GOODS fields, based on simulations aiming to recover art i- 
ficial sources inserted into the images (|Papovich et al.ll2004l ). 

One potential problem with the above matching proce- 
dure is that there could be some instances where a given 
IRAC source could be matched to several GNS galaxies. 
The deblending procedure we use accounts for this problem 
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in most cases. However, we have identified seven instances 
where one infrared detection could be associated with mul- 
tiple HST/NICMOS sources. We have removed these seven 
objects from this analysis so that no ambiguity remains. This 
means that all galaxies in the mass complete GNS sample 
of 1282 galaxies have IRAC counterparts that are secure. 
MIPS detections are not as high, with -21% (267/1282) of 
the GNS galaxies detected above a flux limit of f24 = 30 fj, Jy. 

The total infrared luminosity, Lir (8-1000/im), 
is determined follo wi ng t he pro cedure described in 
IPerez-Gonzalez et alj f2008) and IPerez-Gonzalez et al.l 
( 200&t ) . Using the o bserved 24 p,m f lux d ensity, we use 
the SED templates of lCharv fc Elbaj (|200lh to best repro- 
duce the corresponding rest-frame luminosity at ~ 8/im at 
z — 1.5 — 3, and then obtain the integrated Lir of that 
template. If a source is also detected at 70 /im, then a 
colour is used in the comparison with templates, instead 
of a monochromatic luminosity. We use the Lir in Sec- 
tion 13.31 to determine infrared SFRs. Rest-frame luminosi- 
ties at 1600A and 2800A were also estimated from the in- 
terpolation of the best-fitting templates at these respective 
wavelengths. A UV-slope was then calculated from the inter- 
polated monochromatic luminosities, which will be further 
discussed in Section T3. 2 1 



3 STAR FORMATION ACTIVITY 

A major goal of this paper is to determine, for the first time, 
the star formation activity present in a stellar-mass-selected 
sample of galaxies in the redshift regime 1.5 < z < 3. 
Determining star formation activity at these redshifts is not 
trivial. Infrared observations are useful indicators of dust 
heating due to star formation, but current Spitzer Space 
Telescope observations are not yet deep enough to detect 
the full mass-selected sample of galaxies, as only 60% of the 
massive galaxies are detected, and only —21% (267/1282) of 
the entire sample are detected at 24 /im (above a flux limit 
of f24 = 30 /iJy; see Section [2~5)l . 

On the other hand, 98.5% of the mass-selected GNS 
sample are detected in the rest-frame ultraviolet (observed 
ACS/zgso-band), which is sensitive to young (< 1 Gyr 
old) stars ev e n for high redshift, optically red galaxies. 
iKavirai et alJ |2008) conclude that in their sample of early- 
type galaxies at 0.5 < z < 1, the UV flux comes exclusively 
from young stars and the same holds for our even younger 
galaxy sample. Furthermore, they find that only —1.1 per- 
cent of those galaxies are consistent with purely passive age- 
ing since 2 = 2, providing confidence that the rest-frame UV 
flux from galaxies in this study is highly sensitive to young 
stars, and therefore representative of recent star formation. 
In this section, we derive and compare star formation rates 
from both rest-frame ultraviolet and infrared wavelength 
regimes and determine the best star formation measure to 
use for the rest of this study. 



3.1 Ultraviolet Star Formation Rates 

We use rest-frame UV light as a SFR indicator, which we 
measure from the observed optical ACS/zsso-band. The UV 
luminosity is closely related to the level of ongoing star 
formation as the UV continuum luminosity is produced 
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Figure 2. The difference between A2800,SBD based on an SED- 
determined UV slope and A280O06S1 derived from the observed 
(B — z) colour. 



by short-lived O an d B stars, even at redshifts of z > 1 
l|Kavirai et al.|[2 008). We determine the SFRfy v from the ob- 
served optical ACS/z850-band flux density (with a 5 a limit 
of 27.5 in the AB system), which corresponds to the rest- 
frame UV luminosity around 2800 A, spanning wavelengths 
of 2125 - 3400A for z = 1.5 - 3 galaxies. Of the 1282 GNS 
galaxies in our sample with log Af* ^9.5 and 1.5 < z < 3, 
only 1.5% do not have counterparts in the ACS/zsso-band 
catalog (see Section [2.11) . 

To calculate SFRyv, first we derive absolute magni- 
tudes using the kcorrec t package (Blanton fc Roweis!l2007l . 
v4.2). We then use the lKennicuttl ~( 19981 ) conversion from 



2800A luminosity to SFR assuming a Salpeter IMF: 

SFRuv (M yr" 1 ) = 1.4 x KT 28 L2800 (ergs s" 1 Hz" 1 ). (1) 

Before dust extinction is taken into account, we find at 
z = 1.5 a limiting SFR,yv, i> s = 0.28 ± 0.1 M Q yr _1 and at 
z = 3.0, we find a limit of SFRr/y, & s = 0.98±0.3 M Q yr _1 . 
The errors quoted here take into account photometric er- 
rors and the error in the conversion from a luminosity. The 
overall errors for individual SFRs are about 30%, which are 
dominated by the dust correction applied, as we discuss in 
the next section. 



3.2 Dust Absorption Determination 

As is now well known, SFR estimates derived from UV light 
can be highly affected by dust extinction. Inferring the level 
of dust reddening from photometric properties remains a 
major challenge for understanding the star-forming proper- 
ties of distant galaxies. To address this issue, we examine 
techniques to recover UV light absorbed by dust and re- 
radiated at longer wavelengths. 

One method for correcting for dust obscuration is to 
use t he UV slope a s an e stimate for the amount of extinc- 
tion. iMeurer et alj (Il999l ) demonstrated a correlation be- 
tween dust attenuation and rest-frame UV slope, /3, for a 
sample of nearby starburst galaxies (where Fx — A' 3 ). Up- 
dated studies of local galaxies using the Galaxy Evolution 
Explorer (GALEX) near-ultraviolet ban d dBuat et al.ll2005l ; 
ISeibert et ai1l2005h and z - 2 galaxies (|Reddv et alJ|201Ch 
show that the UV slope from the local starburst relation can 
be used to recover the dust attenuation of a vast majority of 
moderately luminous galaxies at z — 2, b ut with a large scat- 
ter of a — 0.9 mag l|Seibert et alJ koOS). Dust attenuation 
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log SFRuv,corr(obs) [M SU n/yr] redshift 

Figure 3. The left panel compares SFRyy- derived from observed zgso-band with different dust corrections applied. Along the X-axis of 
the left-hand panel is SFR^y jCorr ( f, s ) , corrected by observed optical colours; along the Y-axis is SFRjyy- , C orr(SED)> corrected from an 
SED-determined UV slope. The solid line shows the one-to-one relation while the dotted lines shows factors of 10 difference. The right 
panel shows SFKijv,corr(SED) as a function of redshift. Open circles identify massive galaxies with log Mr ^ 11. Small grey points show 
all GNS galaxies with log M* > 9.5 and 1.5 < z < 3 not detected at 24 fim. Black points show galaxies detected at 24 /im with /24 > 30 
/ijy. Grey arrows show the redshifts of objects fainter than the 5 a zgso-band detection limit of 27.5 mag (AB). Points in the left plot 
at SFKijv,corr(obs) = —0.2 show galaxies below the 5 a _B-band limit. 



measurements derived from the UV for galaxies at z ~ 2 
have been examined for samples of galaxies s elected from 
several different methods : the BzK method ()Daddi et al.l 
120041: IPannel a et al.ll2009l l: Lyman Break Galaxies (LBGs ) 
jReddv et al.ll20ld ); and rest-frame UV l|Erb et al.l 1200^ . 
However, as is well known for red galaxies, sorting out the 
difference between dusty galaxies and those with evolved 
stellar populations can be difficult (e.g. Ivan Dokkum et "all 
l20ld : IViero et al.ll20lJ ). 

Another method for determining dust extinction re- 
lies on using the reddening parameter value extracted from 
a best-fitting SED template. As described in Sections 12.31 
and 12.51 we fit sets of templates built from the stellar popu- 
lation synthesis mo dels twice for this sample: once to derive 
the stellar masses jGriitzbauch et al.| [2010: Conselice et alj 
l201ll) and separately to de termine the infrared lumniosity 
ijPerez-Gonzalez et al.ll200&t ). 

Unfortunately, neither of our SED fits provide reliable 
values of dust extinction and we therefore rely on the UV 
slope method to determine a dust attenuation, A2800- There 
are some notable limitations and caveats that exist when 
using this method to correct for dust, since the technique 
assumes that the UV slope is due to dust reddening, rather 
than other sources such as evolved stellar populations. De- 
spite the difficulties, we are encouraged by the facts that a 
large majority of the GNS sample shows colours of galaxies 
experiencing some level of star formation as in Figure [T] and 
98.5% are detected in the rest-frame UV. 

Here we apply two methods for determining a UV dust 
attenuation, A2800, in terms of the UV slope, as well as 
investigate proper t ies of our red distant galaxies (see also 
iGrutzbauch et alJ l|2010h for a discussion of this). But con- 
sidering the difficulties in determining dust extinctions in 



the UV, for the results presented in this paper, we show 
SFR;yy- before and after correcting for dust. 

The first method we investigate is an empirical tech- 
nique for sta r- forming ga l axies that comes from observed 
optical data. iDaddi et al.l |2004l ') calibrated the correlation 
between E(B — V) and the observed optical (B — z) colour 
using a BzK-selected sample of star-forming galaxies with 
K < 20 in the redshift r ange 1.4 < z < 2.5. Combining this 
with the lCalzetti et all (|2000h reddening law, this becomes: 



= 1.815 {B -2 + 0.1) AB . 



(2) 



We note that this method is only applicable to the 87.5% 
of the mass-selected GNS galaxies that are detected above 
the 5<r Bab limit of 28.2. Applying this method we find an 
average value of A2800, obs of 1.9 ± 1.2 mag. 

The second method we use to determine A2800 utilises 
the SED-fitting procedure described in Section 13.11 Since 
we do not recover a reliable Ay as a free parameter of the 
fit, we use the best-fitting SED to extrapolate and measure 
UV luminosities at wavelengths of 1600A and 2800A, and 
use these to calculate the UV slope /3. This should be a 
more reliable approach to uncovering the dust extinction as 
it includes older evolved SEDs, and provides a value of A2800 
for every galaxy. 

We follow the ICalzetti et "all l^OPd ) methodology 
to derive Aisoo = 2.31/3 + 4.85 and A2soo,sbd as 
A 2 800,sbd =Ai6oo(fc280o/fci60o)- The slope, P, is calibrated 
for the values -2.0 < ^ < 1.0 with the lCalzetti et al.l (120001 ) 
law and for this reason, if the SED-derived /3 is less than 
— 2.0 (144 objects), we assign a value of A2800 = 0.0. If 
/3 > 1.0 (60 objects), we assign an A2800 = 3.5 so that there 
is not a large overcorrection for dust. Both of these cases are 
however rare. A similar upper limit of A.fuv ~ 3 was em- 
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ployed bv lSalim et all (|2007l ) based on the (FUV - NUV) 
colour observed by GALEX. Using this metho d we find an 
averag e value of A2soo,sbd = 1.6 ± 1.2 mag. iTresse et al.l 
(2007) study galaxies from the VIMOS VLT Deep Survey 
(VVDS) and find similar values for far-UV dust obscuration 
of 1.8 — 2 mag from z = 2 to z — 0.4. 

Figure [2] compares the two separate determinations of 
A2S00 as a function of stellar mass and redshift. They are in 
relatively good agreement for M* < IO^Mq with the abso- 
lute value of < AA >= 0.3 and standard deviation of 0.6. 
The scatter increases for M* ^ lO n M0 galaxies showing 
that there is a difference between the empirical and SED 
approach towards dust corrections. Values of A2800, obs tend 
to be larger than A2soo,S-BD for M* ^ 10 1 Mq galaxies by 
the amount 5A2800 = 0.86 ± 0.7. As shown in the right 
panel of Figure [2] the two measures of A2800 are, however, 
in rough average agreement across the entire redshift range, 
with < A A > = 0.18±0.03 for 1.5 < z < 2.6. At z > 2.6, 
the A2800, obs begins to leave the redshift range for which 
it was calibrated and we see a systematic offset such that 
A2800, obs is larger than A2800,s_e d by a value o f 0.8. T hese 
results confirm recent findings bv lNordon et alj (|2010h who 
stack Her schel observa t ions f rom the PACS instrument to 
show that iDaddi et al.l l|2004l ) calibration overestimates the 
dust attenuation. 

In Figure [3] we show how these separate dust correc- 
tions affect the actual SFR measures. The Y-axis shows 
the SFKiJv,corr(obs) corrected with the empirical method 
(A 2 800,of> 3 ) and along the X-axis we show SFR UVtCorr ( S ED), 
corrected using the SED fitting based dust correction 
(A2800,S-E.d). Nearly all measurements of these two dust cor- 
rections agree within a factor of 10 (grey dotted lines). 

There is evidence for an offset present among massive 
galaxies not detected at 24 /j,m (open grey circles with dots) 
in Figure O This difference is such that the dust correction 
derived from the observed UV slope, and thus the resulting 
star formation rate, is larger than the SFR corrected with 
the SED-derived UV slope. This difference can be quite a bit 
larger, over a factor of ten. This demonstrates that the SED 
method for finding the dust extinction, which we ultimately 
use, is able to locate those systems which are not dusty, but 
red due to older stellar populations. The star formation rates 
corrected with the SED-derived dust extinction are also on 
average lower than the empirically corrected SFRs, and are 
thus if anything underestimating the star formation rates 
for these galaxies. 

What is clear from this comparison is that there are red 
galaxies, particularly with large masses, that have older stel- 
lar populations, creating their red colours. The A2800 values 
calculated from the obs e rved (B — z) colour was originally 
derived by IDaddi et af] (|2004l i from a set of galaxies cho- 
sen to be star-forming, and therefore will not account for 
galaxies with old stellar populations. Since the determina- 
tion of the UV slope from observed passbands was originally 
derived from star-forming galaxies only, we apply the SED- 
derived dust correction to our observed SFR[/y,of>s and use 
this as the main SFR indicator throughout the rest of the 
paper. 

Furthermore, the validity of our dust corrections are 
backed up by comparing the fraction of galaxies detected at 
24/^m (§ 13. 3|) as a function of our completely independently 
derived values of /3. The left panel of Figure [4] shows the 
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Figure 4. Fraction of galaxies detected above the limiting 24 (im 
flux of 30 fijj as a function of UV slope 0, stellar mass, and rest- 
frame (U — B) colour. Fractions derived from the full redshift 
range are shown as solid diamonds, while galaxies restricted to 
the lower redshift range of 1.5 < z < 2 are shown as open circles. 



fraction of galaxies detected at 24 ^m as a function of /3. We 
find a step decline in the fraction of galaxies detected from 
~ 60% at P ~ 1.5 down to 15% at P ~ -1.5. This shows 
that our measures of dust extinction, through the value of 
/3 derived from SEDs, correlates well with the presence of 
24/mi emission from heated dust, as would be expected. 

The dust corrected SFKuv,corr that we will use 
throughout the rest of this paper (the rest-frame UV SFR 
derived from the observed ACS/zsso-band, corrected using 
the SED-determined A2800 values) is shown in the right 
panel of Figure [3] as a function of redshift. We find that our 
galaxies exhibit a range in SFRfjy jCorr between 0.6 M0yr _1 
and ~2000 M0yr _1 . The average SF~Ruv,corr for the full 
mass-selected sample between 1.5 < z < 3 is 42^29 M0yr -1 . 
For massive galaxies with Af* ^ W 11 Mq, we find an average 
SFRuv.corr of 103 ± 75 Moyr -1 (consistent across the full 
redshift range). Among galaxies detected at 24 /im (~21% of 
the sample), we find average SFKuv,corr = 73 ±50 M0yr~ . 
The errors quoted are the dispersion of the values of SFR in 
the sample. The majority of SFRs of individual galaxies have 
errors of 0.2 dex, which are dominated by the dust correc- 
tion applied. Furthermore, our value s of A2800 and SFRjjy 
are consistent with those derived by lOnodera etaD (|2010h 
who derive SFR^y from rest-frame 1500A luminosity. They 
observe a sample of star-forming BzK galaxies at z ~ 2 and 
determine dust corrections from optical emission lines and 
test their agreement with the rest-frame UV slope. Over- 
all our dust-corrected values of SFR(yy within their stellar 
mass range of M* > 5 x 1O 1O M0 are consistent with their 
findings. 



3.3 Infrared Star Formation Rates 

Star formation rate and dust properties of distant galaxies 
can also be measured through emission from the infrared. 
This light originates from dust in galaxies which has been 
heated by photons from massive stars, and thus in principle 
should agree with measurements of the star formation rate 
derived directly from these stars in UV light once corrected 
for dust extinction. 

Over the full redshift range from 1.5 < z < 3, we find 
that ~21% (267/1282) of our sample is detected at 24 /im, 
above a flux limit of f24 = 30 /jJy. If we restrict the redshift 
range to 1.5 < 2 < 2 then we find a higher fraction of 35% 
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detected. In Figure [4] we show detection fractions at 24 /j,m 
for galaxies over the full redshift range (solid diamonds) and 
for a subsample of galaxies restricted to a lower redshift 
range of 1.5 < z < 2 (open circles) separately. The trends 
are qualitatively similar for both groups and we give values 
for the full redshift range unless otherwise indicated. We 
find that the detection fraction at 24 /im depends on stellar 
mass and rest-frame (U — B) colour. Nearly 60% of the GNS 
galaxies with log M* > 10.8 are detected at 24 /im. The 
detection fraction decreases to ~20% at log M* ~ 10.25 and 
less than 5% for log M* < 10 galaxies. Between 50 and 60% 
of galaxies with (U — B) ^ 0.8 are detected at 24 )im while 
the detection fraction decreases at bluer colours. A detection 
fraction of 30% is found for galaxies with 0.6 < (U—B) < 0.8 
and 10% for the bluest galaxies at (U - B) < 0.6. This 
implies that more massive galaxies and galaxies with red 
colours have higher obscured SFRs and more dust, causing 
enough emission at 24 fim to be detected. 

We can get some idea about the relative SFRs mea- 
sured from UV and FIR light through Figure [5] which 
presents the ratio of Lir to L2800, as a function of the SED- 
derived UV slope, f3. Galaxies with lower infrared luminosi- 
ties, 11 < log Lir < 11.5, most closel y mimic the local rela - 
tion for starburst galaxies derived by iMeurer et all (QUI). 
As infrared luminosity increases, the Lir to L2800 ratio in- 
creases such that nearly a l l ULI RGs (with log Lir ^ 12) 
lie above the lMeurer et all (I1999T ) relation (right side of Fig- 
ure [5]). 

In our sample of galaxies, we find that 44% of all galax- 
ies detected at 24 /jm are classified as ULIRGs. These re- 
sults largel y agree with the z ~ 2 star-forming galaxies 
studied in [Reddv et all (|201Ch who select UV-bright Ly- 
man break galaxies . Our values have a larger spread than 
I Reddv et all ((2010), most likely due to the fact that we 
are utilising a mass-selected sample. This larger variation of 
star-forming propertie s may be especially promi nent among 
massive galaxies (e.g. Ivan Dokkum etafl l201Ch . Similarly, 
using the GOALS sample iHowell et all (|201Ch find that 
LIRGs have larger IR excesses than lower luminosity galax- 
ies of similar UV colour. On average, we also find that LIRGs 
and ULIRGs have larger IR excesses. 

To de r ive in frared SFRs from from Lir, we use the 
iKennicuttl (| 19981 ) conversion, 



SFR/ij (M Q yr" 1 ) = 4.5 x 10~ 44 Lir (ergs s~ 



(3) 



We then determine the total SFR in the infrared for 
each galaxy by adding the SFR/j? together with the ob- 
served value of SFRuv,uncorr ■ Figure [H] directly compares 
SFR/K+t/v and SFRuv.corr for mass-selected GNS galaxies. 
Among massive galaxies (open circles), ~87% agree in SFR 
within a factor of 10, while 80% of all the galaxies in Fig- 
ure|6]agree to within a factor of 10. However, there is no real 
correlation between these measures and the scatter is large. 

We find that the SFR/ij+yv is an average factor of 
3.5 Ij'j larger than SFRuv.coir between 1.5 < z < 3, an 
offset whi ch has been seen befo r e in galaxie s at h igh red- 
shift (e.g. iPapovich et al. 1 120071 ; lElbaz et all |2010| ). There 
are several fact ors that could con tribute to this offset. A 
recent study bv lElbaz et ail (|2010h compares Herschel deep 
extragalactic surveys from 100 to 500 )im total IR emission 
extrapolated from 24 /im Spitzer fluxes. They estimate that 
at z > 1.5, the 24 /im observations overestimate the SFR by 
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Figure 6. SFRjfl + [/v' versus SFRuy corr for all galaxies de- 
tected at 24 fim. Open circles identify massive galaxies with 
M, > 1O 11 M0. The solid line shows the one-to-one relation while 
the dotted lines shows factors of 10 difference. 



a factor of ~3, which could be due to a larger PAH emission 
in distant galaxies, hot dust heated by a buried AGN, or 
uncertainty in local SED templates. 

Some of these issues are explored in a forthcoming study 
bv lWeinzirl et all (|201lT ) of massive GNS galaxies with M» > 
5 x 10 10 Mq in the redshift range 1 < z < 3. They determine 
SFR/_r for their sample independently of this study using 
infrared luminosities derived from 24 /j,m data and apply 
the factor of ~3 correction sugg ested by the early result s 



from Herschel (e.g. lElbaz et al]|2010l ; iNordon et all l2010h 
For massive galaxies, they find that the SFR/_r agrees with 
the SFR(7i/,con- values within a factor of ~ 2 to 3 for galaxies 
with SFR/ii > 40Moyr -1 , but galaxies with low SFR/ij are 
lower than SFH,uv,corr by up to a factor of 10. 

It is also possible that the UV-based SFRs could be un- 
derestimated in cases of very high extinction, where star for- 
mation could be completely unrecoverable by UV observa- 
tions alone. It is unlikely, however, that we are missing many 
star-forming galaxies altogether from our samples though, 
since the primary selection is done from deep H-haxid imag- 
ing and a high level of extinction would be needed for an 
object to be comp letely undetectable (see Section [2.41 and 
iFraver et all 120041 ). We are reassured by the range in val- 
ues we find for the r atio between these SFRs are consistent 
with those found in iDaddi et all (I2007I ). who studied BzK- 
selected galaxies between 1.4 < z < 2.5 in the GOODS 
fields, and who also compared infrared and UV star forma- 
tion rate indicators. 

Since the main conclusions we draw in this work are 
based on the mass-complete GNS sample using UV-based 
SFRs, and the corrections derived to apply to the 24 /im 
data are still new, we do not attempt to make any correction 
to the IR-based SFRs in this study. 
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Figure 5. The ratio of L /J? (8-1000Atm) to L f/y (2800A) as a fun ction of the UV-slope, j3, for 1.5 < z < 3.0 galaxies with f24 > 30/iJy. 
The solid line shows the local relation from iMeurer et al ] {Hal). Black points in the right panel show galaxies classified as ULIRGs 
(log Lib > 12). In the left panel, galaxies with 11.5 < log L/u < 12 are shown as red diamonds and galaxies with 11 < log Ljjj < 11.5 
are shown as blue squares. 



4 RESULTS 

In this section, we discuss the star-forming properties of 
GNS galaxies as a function of stellar mass and rest-frame 
(U — B) colour. We discuss the effects of the dust determina- 
tion and the photometric redshift errors on these relations. 
We also investigate what properties may account for the red 
colours observed in galaxies at high redshift. 



4.1 Star Formation and Stellar Mass 

Several studies in recent years have investigated the change 
of SFR as a f unction of stellar mass and time to redshifts of 
z ~ 1.5 (e.g . iBrinchmann fc Ellis! 1200(3; iBauer et all 120051 ; 
iBundv et all 120061 ; IZheng et all 120071 ; INoeske et all 120071 ) 
and have found a general trend for SFRs to rise with increas- 
ing stellar mass. Studies of galaxies at higher redshift, be- 
tween 1.5 < z < 2.5, using SFR indicators from many differ- 
ent wavelength regions are also beginnin g to show evidence 
for increasing SFRs with stellar mass (e.g.|Paddi et al.ll2007l ; 
lRenzinill2009l ; ISantini et alll2009l ; lOnodera et"alll2010l ). We 
present here for the first time the SFR-stellar mass relation 
for a sample complete in stellar mass to M t = 3 x 1O 9 M 
between 1.5 < z < 3. 

Figure [7] shows SFRuv,corr as a function of stellar mass 
in three redshift bins, from z — 1.5 to z = 3. It is clear 
that a trend between SFR and stellar mass exists at each of 
these redshifts. Within a single redshift bin, the median SFR 
increases by roughly a factor of 10 from the lowest stellar 
mass galaxies probed, log Af, = 9.5, to the highest mass, 
log M» ~ 12 galaxies. Galaxies not detected at rest-frame 
UV wavelengths are considered non-star-forming. Non-SF 
galaxies compose only 1.5% of the sample and are shown as 
downward facing grey arrows in Figure [7] In computing the 
median SFR[/y iCorr values (large blue circles), we include 
galaxies not detected at rest-frame UV wavelengths in each 
mass and redshift bin. 

The correlation between stellar mass and SFR for the 



GNS sample has a slope consistent with that found at inter- 
mediate redshift bv lNoeske et al. (2007), shown as the green 
dashed line in Figur e [7] Using the same stellar mass range as 
INoeske et "all l|2007l ) (10 < log M, < 11), the normalisation 
found for GNS galaxies at z ~ 2 is a factor of 2.2 larger than 
at z ~ 0.7. At fixed stellar mass, star-forming galaxies are 
more actively forming stars, on average, at z ~ 1.5 than at 
z ~ 0.7, lik ely due to a larger a bundance of gas that depletes 
with time (|Noeske et al.ll2007l ). 

In a study of if -selected GOODS galaxies at z ~ 2, 
iDaddi et all i|2007l ) also find increasing SFRs with stellar 
mass, as shown by the dash-dot line in the top row of Fig- 
ure [7] Th e correlation of the GNS sample has a flatter slope 
than the IDaddi et all i|2007h relation at z ~ 2, but simi- 
lar median SFR values for log M* < 11 galaxies. Among 
massive log M, ^ 11 GNS galaxies, we find an average 
SFKuv,corr = 103 ± 75 Mo yr" 1 , which is roughly a fac- 
tor of 1.6 less than the lDaddi et all (120071) relation . We find 
a flattened relation relative to Daddi et all |2007l ) which is 
either due to the fact that we are using a mass-complete 
sample instead of just star-forming galaxies, o r due to the 
overes timation of the dust correctio n appl ied bv lDaddi et all 
(2007), as found by iNordon et all (2010) and discussed in 
Section |3 . 2 1 ab ove . Interestingly, for a given stellar mass, we 
see roughly constant SFRs for all galaxies over the redshift 
range 1.5 < 2 < 3 for the full mass-selected GNS sample. 

We plot the 19 non-star-forming galaxies as grey arrows 
in Figure [7] Recall from Section 13.11 that non-star-forming 
galaxies are those not detected at the observed ACS/zsso- 
band 5 a limit of 27.5 in the AB system, corresponding 
to an uncorrected SFRyy of 0.98 ± 0.3 Moyr -1 at z = 3. 
Only one non-star-forming galaxy is at 2 < 2 and only four 
are detected at 24 /im, implying that we are not missing a 
significant population of highly extincted starbursts at z > 
2. The non star-forming galaxies are mostly low mass, and 
only three have M, > W^Mq. We will look more at the 
colours of these galaxies later, in Section |4. 41 

We present the specific star formation rate (SSFR = 




log M/Msun 

Figure 7. Median SFRy v,corr versus stellar mass in three redshift bins (top) and median SSFR versus stellar mass (bottom). Median 
SFRs are shown for GNS galaxies in each stellar mass bin as large blue circles. Individual GNS galaxies are black points, with galaxies 
below the 2:g5o-band detection limit shown as grey arrows. Grey diamo nds show results fro m lDrorv fc Alvarej l|2008h who use SED fitting 
to derive UV SFRs. The dash-dot diagonal line in t he top row is thelDaddi et al. relation for galaxies at z ~ 2 while the green 

dashed line shows the intermediate relation found bv lNoeske et al] d2007T> at z ~ 0.8. 



SFR/M.) as a function of stellar mass in the bottom three 
panels of Figure [7] The error bars show the standard devia- 
tion of the distribution of galaxies in each bin. The median 
SSFR of galaxies with log M* ^ 10.5 is largely independent 
of stellar mass and redshift over the ~ 2 Gyr period between 
z = 3 and z = 1.5. Galaxies with log M* 10.5 show a 
sharp decrease in SSFR with increasing stellar mass such 
that the median SSFRs decrease by 0.6 to 1.0 dex between 
log M» ~ 10.5 and log M» ~ 11.5. 

Using radio observations of SSFRs, iPannella et al] 

(2009) also find that the average SSFR is constant 
over almost on e dex in st el lar m ass at z = 1.6 
and z = 2.1. iDunne et al] ((2009) also find a shal- 
low decrease of SSFRs with decreasing stellar mass us- 
ing stacked radio observations. As a function of redshift, 
we see no dependence of SSF R with redshift, in con- 
trast to studies at lower redshift llBrinchmann fc Ellisll2000l; 
Bauer et al.ll2005l; iPerez-Gonzalez et al] 120051 ; IZheng et al] 
20071 ; I Pannella et alj|2009l ). but in agreement with recent 
results from the infrared using the Hersh el Space Telescope 
l|Rodighiero et al]|201fj ICava et aljboich . 

It is possible that the low SSFRs we see at z ~ 3 for high 
mass galaxies is indicative of a period of stellar mass growth 
dominated by merging events, which subsequently slow at 
z < 3 so that star formation becomes the dominant growth 



mechanism for the most massive galaxies until their cold gas 
supply runs out. A massive galaxy shown in Figure [7] with a 
SSFR between the median values of logSSFR = —9.0 yr _1 
and —9.5 yr" 1 will double its stellar mass from star forma- 
tion alone in 1 to 3 Gyr. 



4.2 The Effect of Photometric Redshifts 

The errors on the photometric redshifts determined for GNS 
galaxies are consistent with those found in the literature for 
galaxies in this redshift range for studies using a similar 
number of filter bands, as described in Section 12.21 but in 
order to further test the effect of the photometric redshift er- 
rors on our results, we have perf ormed a set of Monte Car lo 
simulations (fully described in iGriitzbauch et all l|2010h ). 
We randomise the photometric redshift input according to 
the typical photometric redshift error in our redshift range 
of Az/(1 + z) ~ 0.1, obtained by the comparison to the 
available secure spectroscopic redshifts for the galaxies in 
our sample. 

We assume a Gaussian distribution of errors on the 
photo-z redshift for each object, where the width of the dis- 
tribution corresponds to Az/(1 + z) ~ 0.1 for all galaxies in 
the redshift range 1.5 < z < 3. For each galaxy, a random 
value is selected within this distribution, which is then added 
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Figure 8. Median SFR versus stellar mass in three redshift bins (top) and median SSFR versus stellar mass (bottom). The blue points 
are the same SFHjjv,corr as in Figure [7] Light blue inverted triangles show SFRuv, uncorr for all GNS galaxies. Open red squares are 
median SFRj^+jy^/ for GNS galaxies detected at 24/^m. Open gray triangles show SFRr^; from Perez-Gonzalez et al. (2008). 



5^ 

c 



oc 

LL 

CO 
CO 

o 



I 

DC 

LL 

CO 
CO 

CD 
O 



■ | I | I | I | I | I | I | I | I | I | I | 

1.50<z<2.00 

® 

©. ® 



-10 - 



w; . « « 
.*. ■„ " • • 



I I' 1 ' 

1 .50<z<2.00 , 

- :• '- : '®: .. ~r 





ULIRG ©©•• 
o log M/Msun > 1 1 -6 



■ 1 ■ 1 ■ * ■ 1 ■ 1 ■ 1 ■ 1 ■ * ■ 1 ■ 1 ■ 1 ■ 



■ I i I i | i I i I i I i I i | i 

2.00<z<2.50 







"79:1, •■ 



I 

2.00<Z-<2.50 

x 

• ■ - 





• ©^ 

■ 1 ■ 1 ■ I ■ 1 ■ 1 1 1 ■ 1 ■ I ■ 



I I I I I I I I I I I I I I I I I I I I I 

2.50<z<3.00 

® 

o 

.... ® 

® I 

>H>© 



Mi 



© 



I I' 1 

2.50<:z<3.00 

' 4PP* 

© 

■ 1 1 1 ■ I ■ 1 ■ 1 ■ 1 ■ 1 ■ I ■ 1 ■ 1 ■ 1 ■ 



0.5 1.0 



0.5 1.0 



0.5 1.0 



(U-B)restframe 

Figure 9. Star formation rate versus rest-frame (U — B) colour in three redshift bins (top) and specific star formation rate versus 
(U — B) colour (bottom). Red crosses show galaxies classified as ULIRGs (L/^ ^ 10 12 Lq) and open circles identify massive galaxies 
with M» > 1O 11 M . 
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to the measured photo-z. Galaxies outside our redshift range 
of interest are included in the randomisation to account for 
scattering in and out of our redshift range. Catastrophic 
outliers are accounted for by randomly adding much larger 
offsets to the original redshifts, of the percentage of galaxies 
corresponding to the fraction of catastrophic outliers, which 
is calculated from the comparison with secure spectroscopic 
redshifts. 

We then repeat the computation of the SFR and SSFR 
as a function of stellar mass and colour with the new photo- 
metric input. The results indicate that while the individual 
values of SFR and SSFR for a single galaxy can change sig- 
nificantly, on average, the relations we report in the Sec- 
tion 14. f I and later in Section 14.41 remain robust and un- 
changed despite the errors of the photometric redshifts used. 

In order to determine how much the lower photomet- 
ric redshift reliabilities of fainter galaxies affect the results 
presented here, we repeat the Monte Carlo simulations us- 
ing a magnitude-dependent approach. We calculate the pho- 
tometric redshift error and outlier fraction separately for 
the bright and faint galaxies separating them at H = 24 
mag. The main impact of this procedure is an increase in 
the catastrophic outlier fraction for the faint galaxies in our 
redshift range to 33%, as opposed to 17% for the brighter 
galaxies. The random error stays about constant for both 
samples at Az/(l + z) ~ 0.1, whereas the offset increases 
slightly from ~0.05 for the brighter sample to ~0.08 for the 
fainter galaxies. 

The result of repeating the simulations with the above 
values separately for bright and faint galaxies is that the 
number of galaxies in each bin decreases due to the large 
outlier fraction which scatters one third of the faint galaxies 
out of the redshift range of this study entirely, but we find 
that the overall results remain unchanged. We do not correct 
the measured values as a result of the Monte Carlo simula- 
tions, because we have found that the uncertainty around 
the resulting relations is hardly affected, other than the scat- 
ter increases slightly. While the uncertainty is larger for less 
massive galaxies, the results are robust to various photomet- 
ric redshift errors and improved reliability will be a goal of 
future spectroscopic redshifts, but that these will have to be 
done with 20-30m telescopes due to the faintness of these 
systems. 

4.3 The Effect of Dust 

In order to see how the influence of dust affects these results, 
we show several different measures of SFR in Figure[H] which 
shows the SFR and SSFR as a function of stellar mass. The 
blue circles are the same SFRy v.corr presented in Figure 
Light blue inverted triangles show the SFRpv,„„ cor r before a 
dust correction is applied. The SFRj7y, uncor . r are consistent 
with a slight increase with stellar mass in each redshift bin, 
and also a slight increase in SFR as a function of redshift in 
each stellar mass bin, except for the most massive galaxies 
which have the same median values at all redshifts. This 
implies that the dust cor rection is a function of stellar mass 
(|Criitzbauch et aHl2010Ti and not only steepens the SFR-M, 
relation but also flattens the SSFR-M* relation. 

The open red squares in Figure [5] show the median 
SFR/_r+!7v for GNS galaxies detected at 24 /im. At z < 2.5, 
the SFR/ij+(7v and SFKuv.corr agree quite well. On the 



other hand, at z > 2.5 the SFR/ij+yy is greater than 
SFKuv.corr by an average factor of 5 and greater than 
SFKuv,uncorr by an average factor of 5 0. A high redshift 
infrared study by IPerez-Gonzalez et al.l (|2008h uses Spitzer 
Space Telescope observations at 24fim to determine SFRs 
of galaxies at z ~ 2. Their results for log M* ^ 10.5 
galaxies are shown in the bottom panels of Figure [8] as grey 
triang les connected by a dash-do t line. At z < 2, the SSFR 
from lPerez-Gonzalez et all (|2008T ) are in excellent agreement 
with both SSFRiH+c/v a nd SSFRc/y, corr of the GNS sam- 
ple. As redshift increases, iPerez- Gonzalez et"afl (|2008h find 
higher values of mid-IR-derived SSFR and by z > 2.5, both 
measurements of mid-IR-derived SSFR are one dex higher 
than the SSFRjyv.corr value and nearly two dex higher than 
SSFRc/y, 

uncorr • 

Based on the overall agreement between the median 
dust-corrected SFR^v and the mid-infrared derived SFR 
at z < 2.5 shown in Figure [H] it appears we are successfully 
able to recover the dust obscured star formation using the 
rest-frame UV slope method to determine the dust extinc- 
tion. The large difference at z > 2.5, however, is likely a 
combination of the SFRi_r+[/v being exceptionally hi gh at 
z > 2.5 (e.g. iPapovich et all [20071 ; felbaz et~aH |2010h . the 
SFKuv.corr representing the full mass-selected sample, in- 
cluding galaxies not detected at 24 /im, and the uncertainty 
in using the UV spectral slope to determine UV dust extinc- 
tion at these high redshifts. 

4.4 Star Formation and Restframe Colour 

Based on optical and near-infrared colours of galaxies, there 
is evidence that massive galaxies with red colours at z ~ 2 
are probable progenitors of red non -star-forming massive 
galaxies found a t low redshift (e.g. IConselice et al. 1 120071 ; 
iKriek et al.ll2008T ). Usually, it is assumed that red massive 
galaxies seen in the distant universe consist of old stellar 
populations which passively evolve until the present day. 
However, our results show that it is possible that these mas- 
sive galaxies are red due to dust, or a combination of older 
stellar populations and dusty star formation, as we investi- 
gate in this section. 

One way to identify galaxies with large amounts of 
dusty star formation is to look at populations with excess 
infrared emission such as ULIRGs, classified as galaxies with 
log (Lir/ Lq) > 12. The infrared continuum of a ULIRG is 
dominated by thermal re-radiation by dust. Dust particles 
absorb UV light, mostly emitted by massive main-sequence 
stars in regions of intense star formation activity, and re- 
radiate the energy at longer wavelengths. We saw in Figure[5] 
that ULIRGs have excess IR luminosity relative to UV lumi- 
nosity, and we investigate here whether these galaxies have 
higher dust-corrected SFRs, and from what their rest-frame 
(U — B) colours might originate. 

In Figure O we present SFKuv.corr (top) and SSFRs 
(bottom) as a function of rest-frame (U—B) colour, split into 
three redshift bins. Symbols are the same as in Figure[7|with 
the addition of red crosses indicating galaxies detected at 
24/im that are classified as ULIRGs. It is clear from Figure[9] 
that low-mass galaxies have bluer colours with the majority 
exhibiting rest-frame (U — B) < 0.8. Massive galaxies are 
almost all red with rest-frame (U — B) > 0.8. 

We find that galaxies classified as ULIRGs span the 
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Figure 10. Reddening factor A28OO as a function of rest-frame (U — B) colour (left), stellar mass (center), and SFRuY.corr (right). The 
average A2800 value for all GNS galaxies is shown by the dash-dot line. The open circles show galaxies with log M* ^ 11 while the red 
crosses identify galaxies classified as ULIRGs. 



entire range of rest-frame (U — B) colours but show a trend 
for exhibiting higher SFRs and SSFRs at z > 2. Galaxies 
classified as ULIRGs are believed to be the m ost common 
mode of star formation in massive galaxies (e.g. lDaddi et al.l 
120071 ). The blue dashed line in Figure[9]shows the 2 a clipped 
least squares fit to all the galaxies in the mass-complete 
GNS sample. In the upper panels showing the SFR, we find 
that 80% of the 50 ULIRGs in the middle redshift bin of 
2 < 2 < 2.5 lie above the least squares fit. The fraction is 
59% in the lowest redshift bin and 70% in the highest. We see 
a general trend in the lower panels of Figure [9] where SSFRs 
decrease with redder colours, as seen previously at 2 < 1.5 
jCooper et al.ll2008l ; ls"alim et al.ll2009l ). All the galaxies with 
the highest SSFRs have blue rest-frame (U — B) colours and 
lower stellar masses. This is evidence for the presence of 
dusty star formation in galaxies of all stellar masses above 
the 24 /im detection limit imposed here. 

4.5 What Causes 2 > 1.5 Galaxies to be Red? 

We found in previous sections that massive galaxies are red- 
der, and form stars at higher rates than lower mass galaxies. 
In order to understand what dominates the cause of the red 
colours of these galaxies, we look more in depth at the red 
galaxy population. We define red galaxies to be those with 
rest-frame (17 — B) ^ 0.8 which roughly divides the sam- 
ple into a blue cloud and a red sequence (see Figure [9] and 
Grutzbauch et al. 2011). From the complete mass-selected 
GNS sample, we find that ~ 14% ± 1% (184/1282) are red 
according to this criterion, and among these red galaxies, the 
vast majority, ~ 88% ± 7% (163/184) have log M, ^ 10.5. 
Since massive galaxies dominate the red population, and be- 
cause the 24/im detection fraction drops drastically to less 
than ~25% for galaxies with log M* < 10.5 (see Sec- 
tion [3]3] and Figure [4]), we focus the results in this section 
on log M* 10.5 galaxies. 

In order to distinguish dusty red galaxies, as opposed to 
galaxies that are red due to evolved stellar populations, we 
use three distinguishing characteristics of the galaxy popu- 



lation: (1) the level of dust extinction based on A2800, (2) 
the infrared properties, and (3) star formation rates. 

First, we examine the dust absorption of UV light 
within red galaxies based on the A2800 values derived from 
the SED-determined rest-frame ultraviolet slope. Using this 
method, we denote dusty red galaxies as those that have 
extinction values greater than the average for all galaxies. 
The average extinction value, A2800 = 1.6, is shown as the 
horizontal dash-dot line in Figure 1101 In this way we find 
that among log M* ^ 10.5 red galaxies, a high fraction, 
~ 87% ± 7% (142/163), are red due to dust, and a similarly 
high fraction for log M* ^ 11 galaxies: 89 ± 7%. 

Galaxies with red colours and A2800 ^ 1.6 are not con- 
sidered dusty based on this prescription and therefore likely 
derive their red colour from old stellar populations. This 
amounts to ~ 13% ± 3% (21/163) of the log M* > 10.5 
galaxies and ~ 11% ± 4% (10/90) among the most massive 
log M* ^ 11 galaxies. In this population of massive red, 
non-dusty galaxies, half (5) are curiously detected at 24 /im 
and two of those are classified as ULIRGs, pos sibly due to 
dust heating by another source such as an AGN (|Bluck et all 

ion]). 

Dust attenuation increases as a function of stellar mass 
as shown in t he middle panel of Fig ure 1101 Using the COS- 
MOS survey, iPannella et alj (|2009l ) publish a relation be- 
tween A1500 and stellar mass. Assuming the Calzetti dust 
law to convert their relation to A2800, we find that the ex- 
tinction values agree at log M* = 10.2, but their rela- 
tion is a steeper function of mass such that their A2800 at 
log M» > 11 is higher than the one found in this work by 
roughly one magnitude. This difference is true for the full 
GNS mass-selected sample or just galaxies restricted to the 
redshift range of 1.5 < 2 < 2. This could be due to the fact 
that in this study we assign extinction values of A2800 = 3.5 
to galaxies with /3 > 1.0 (60 objects)in order to avoid over- 
correcting for dust (see Section T3.2p . 

Next, we identify dusty red galaxies as those detected 
at 24/im and classified as ULIRGs, since detection at 24/im 
implies that a galaxy has a significant amount of dust. In 
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terms of 24^m detections alone, the maximum percentage 
of non-dust-dominated red galaxies with log M* ^ 10.5 cor- 
responds to the number of red massive galaxies (163) not 
detected at 24^m, which is 45% ± 5% (74/163). This per- 
centage is quite different than the value of ~13% found using 
the extinction criterion described above. Placing significance 
on whether a galaxy is classified as a ULIRG is slightly ar- 
bitrary, but at z ~ 2, ULIRGs are proposed to be a common 
mode of star formation in massive galaxies (e.g. lDaddi et al.1 
We find that among the 24 /im detections of massive 
red galaxies, 52% ± 6% (47/89) are classified as ULIRGS in- 
dicating that they harbour significant amounts of dust, and 
reassuringly, 92% of those also have above average values of 
A28QO = 1.6. 

As shown in FigurefJ] higher fractions of galaxies are de- 
tected at 24 /im in the redshift range 1.5 < 2 < 2 since we are 
possibly only detecting the most extreme infrared objects at 
2 > 2. If we do the same exercise examining dust in galax- 
ies based on infrared properties with just galaxies between 
1.5 < z < 2, we find that the maximum percentage of non- 
dust-dominated red galaxies with log M* 10.5 corresponds 
to the number of red massive galaxies (79) not detected at 
24/im is 41% ± 4% (32/79), very similar to the value found 
for the full redshift range. The fraction of ULIRGS among 
the low redshift massive red galaxy population is lower than 
for the full redshift range at 43% ± 3% (20/47). 

Despite the red colours of massive galaxies, the aver- 
age SFHuv,corr of red galaxies with log M* ^ 10.5 is 
78 ± 58 Mq yr _1 , consistent over the entire redshift range 
we probe. Recall that among the full GNS sample, we find 
that 1.5% are not detected in the zsso-band and are there- 
fore considered non-star-forming. Among massive red galax- 
ies, we find that 3.6% ± 1.4% (6/163) are not detected in 
the observed zsso-band. If we look to the bluer wavelength 
of the observed B-band (rest frame UV), we find that 87.5% 
of the total sample is not detected at the 5a level. Among 
massive red galaxies, 42% ± 5.0% (68/163) are not de- 
tected in B-band and yet massive galaxies have the highest 
SFRs in the sample and still have red colours, as shown in 
Figure [§] If we assume that dust abs orption is what ca uses 
the lack of detection in the B-band jDaddi et al. 1 12004 see 
also Section [3.2|) . then we can conclude that between 45% 
and 85% of massive red galaxies harbour evidence of dusty 
star formation. If we were overestimating the derived value 
of dust extinction for any of the galaxies in the sample, the 
effect would be to slightly lower this percentage of massive 
red galaxies with dusty star formation. 

We can compare these results to previous work at lower 
redshifts. Previously, it has been shown that extinction con- 
trib utes to the red co lours of galaxies at z = 1. For exam- 
ple, (SmaiTetjil] (|2002r i study a complete sample of 68 EROs 
with (R — K) 5.3 and K < 20.5. Using radio maps and 
photometric classifications, they find that between 30% and 
60% of the ERO population are dusty, star-forming galax- 
ies, as opposed to passively evolved stellar systems. The K- 
selected K20 Survey finds among the ERO population at 
2 ~ 1 that the proporti on of dusty starbur sts at K < 19.2 
is between 33% - 67% (|Cimatti et al.ll2002l ). Recent results 
from the large area COSMOS survey also indicate that half 
of the ne arly 5000 EROs in their sample are dusty starburst 
galaxies (|Kong et al.l 120091 ) based on three different meth- 
ods of classification: SED fitting, infrared colour, and mor- 



phology. Furthermore, IWhitaker etafl (|2010h study massive 
(log M» ^ 11) galaxies in the NEWFIRM Medium-Band 
Survey to 2 ~ 2.2 and find a large scatter in the colours of 
quiescent galaxies. This scatter in rest-frame (U — V) colours 
of quiescent galaxies arises from the spread in ages, due to 
the presence of both relatively quiescent red, old galaxies 
and quiescent blue, younger galaxies toward higher redshift. 

The conclusion is that while massive galaxies have on 
average higher SFRs than lower mass galaxies, they also 
have dustier star formation regions at z > 1.5, causing their 
overall colours to be redder than lower mass galaxies. One 
possibility is that we may be seeing evidence for massive 
galaxies producing more high mass stars due to having a 
more top-heavy initial mass function (IMF). It has been sug- 
gested that discrepancies between the global star-formation 
history and the stellar mass history can be reconciled by 
using evolving IMFs, such that more top-heavy IMFs are 
more common within star forming galaxies at 2 > 2 (e.g . 
iPerez-Gonzalez et al.|[2008l ; IWilkins et ai1l2008l ; lDavell2008r i . 
The most massive stars, i.e., Wolf-Rayet stars create dust, 
and therefore more dust would be produced in top-heavy 
IMF star-forming regions. Another possibility is simply that 
massive galaxies are able to retain more of the dust from 
high mass star formation due to their large gravitational 
potentials, regardless of the IMF. Although, it is possible 
that if we used a more top-heavy IMF to determine stellar 
masses, some fraction of these galaxies would no longer be 
considered massive (M 4 ^ 10 Mq) in our sample. 



5 CONCLUSIONS 

We present the first results of a study designed to in- 
vestigate the star-forming properties of a stellar mass se- 
lected sample of galaxies in the GOODS NICMOS Survey 
(GNS) with masses log M* 9.5 in the redshift range 
1.5 < 2 < 3.0. The GNS is based on deep Hubble Space Tele- 
scope F160W (NICMOS/H-band) imaging of the GOODS 
North and South fields. We combine the GNS NICMOS 
imaging with HST/ACS and Spitzer GOODS photometry 
to calculate photometric redshifts, stellar masses, rest-frame 
colours, dust extinctions, and star formation rates from the 
ultraviolet and infrared. Our sample includes ~1300 galax- 
ies with log M, ^ 9.5 between 1.5 < 2 < 3, of which eight 
percent are massive galaxies with M* ^ 1O 11 M0. 

We measure both UV and IR based star formation 
rates for these galaxies. We determine SFRyy from the 
observed optical ACS 2gso-band flux density which corre- 
sponds to the rest-frame UV luminosity around 2800A. We 
find that the best way to determine the amount of dust 
extinction uniquely for each galaxy is to calculate the ul- 
traviolet slope from UV luminosities calculated from multi- 
wavelength SED-fitting. The average value of A2800 across 
the entire sample is 1.6 ± 1.2 mag. We find that 21% of the 
sample are detected at 24/im above the 30 fiJy (5a) limit, 
and that the detection fraction is much higher, ~ 60%, for 
log M* 11 galaxies. 

The dust-corrected UV star formation rates 
(SFRuv,corr) for all GNS galaxies cover a range between 
0.6 M yr _1 and ~2000 M yr _1 . The average SFRuv.corr 
for the full mass-selected sample is 42 ± 53 MQyr -1 , and 
for just the massive galaxies (M, ^ 10 n M Q ), we calculate 
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an average SFKuv,corr of 103 ± 75 Mqvi- -1 , consistent 
across the entire redshift range. We compare SFRuv.corr 
with SFR/a+c/v determined from the IR luminosity at 
24/im added to the uncorrected UV derived star formation 
rates. Among massive galaxies, ~87% agree in SFR using 
these two methods to within a factor of 10, while 80% 
of all the galaxies in Figure [6] agree to within a factor of 
10. On average, SFR/ij+yy is a factor of 3.5 larger than 
SFKuv,corr between 1.5 < z < 3, with a scatter of 4.6. 

We find that SFRs increase with stellar mass over the 
full redshift range, but the average SFR at a given stellar 
mass remains constant over this ~ 2 Gyr period, opposite 
to what is seen at < z < 1 when massive galaxies have a 
declining SFR. This implies that whatever effect produces 
the downsizing seen in massive galaxies at z < 1 has not yet, 
or has not yet had enough time, to significantly decrease the 
star formation rate within these systems at 1.5 < z < 3.0. 

Furthermore massive galaxies (log M* ^ 11) in our sam- 
ple have the highest SFRs, yet exhibit red rest-frame (U—B) 
colours at all redshifts between 1.5 < z < 3. We conclude 
that between ~45-85% of massive galaxies at these redshifts 
are undergoing dusty star formation, despite a higher frac- 
tion (>85%) exhibiting red colours. We suggest that the red 
colours of massive galaxies is evidence for these systems hav- 
ing dustier star formation regions than galaxies with lower 
stellar mass. One cause for this could be that high mass 
galaxies may have more top-heavy IMFs than lower mass 
galaxies, thus producing more Wolf Rayet stars which are 
responsible for producing dust. 
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